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Abstract 

We study the finite temperature electroweak phase transition with lattice perturba- 
tion theory and Monte Carlo techniques. Dimensional reduction is used to approximate 
the full four-dimensional SU(2) + a fundamental doublet Higgs theory by an effective 
three-dimensional SU(2) + adjoint Higgs + fundamental Higgs theory with coefficients 
depending on temperature via screening masses and mass counterterms. Fermions con- 
tribute to the effective theory only via the Np and mtop dependence of the coefficients. 
For sufficiently small lattices {N^ < 30^ for mn = 35 GeV) the study of the one-loop 
lattice effective potential shows the existence of the second order phase transition even 
for the small Higgs masses. At the same time, a clear signal of a first order phase 
transition is seen on the lattice simulations with a transition temperature close to 
but less than the value determined from the perturbative calculations. This indicates 
that the dynamics of the first order electroweak phase transition depends strongly on 
non-perturbative effects and is not exclusively related to the so-called (j)^ term in the 
effective potential. 
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1 Introduction 



The interest in the study of the high temperature phase transitions in gauge theories, 
initiated a long time ago in has been revived now in connection with possible 
generation of the baryonic asymmetry of the universe at the electroweak scale (see, e.g., 
reviews p|, |^ and references therein). The common belief, based on the perturbative 
calculations of the effective potential for the scalar field, is that the electroweak phase 
transition for moderate Higgs boson masses is weakly first order (i.e. temperature 
metastability range is small compared with the critical temperature). An incomplete 
list of references containing perturbative analysis of the phase transitions is contained 
in Q-||T^. All perturbative calculations, however, work only for sufficiently large values 
of the scalar field, namely ^ gT. This estimate arises as follows. Roughly speaking, 
there are at least three relevant mass scales at high temperatures. The first one is 
related to the Debye screening of the gauge charge and is of the order of = CdqT 
(for the electroweak theory with Np = 6 C/j = ^). The second one is associated 
with the non-abelian magnetic sector of gauge theories, itim = Cmq'^T {Cm is some 
number yet to be determined) and the third one with the gauge boson mass induced 
by the Higgs mechanism, mw = \g4' with being the condensate of the scalar field. 
The naive loop expansion for the effective potential works well provided my/ S> m^i, 
i.e. for (f) ^ T . This condition is usually not satisfied for the interesting range of 
fields and temperatures. The improved loop expansion which takes into account Debye 
mass effects in all orders of perturbation theory works for a larger interval of the scalar 
condensate. However, due to the infrared problem in the high temperature dynamics 
of gauge theories [1^ any of the loop expansions break down for mw < and 
(f) < 2CMgT. Due to the fact that the coupling constant in electroweak theory is not 
so small ((7 ~ 2/3) it is important to know the numerical value of the coefficient Cm- 
If it is small, then perturbative calculations of the effective potential are valid in an 
interesting range of parameters, while a large value of Cm would imply the lack of any 
knowledge of the dynamics of the phase transition. 

The appearance of the magnetic mass is a purely non-perturbative phenomenon 
and not a lot is known about its value. It is associated with the confinement scale 
of 3-dimensional gauge theory derived by dimensional reduction from 4-dimensional 
theory at high temperatures. Different correlation lengths in 3-dimensional SU(2) 
gauge theory were studied by lattice Monte-Carlo methods in JTBI. It was found there 



that the 3-dimensional 0"^"*" glueball mass (inverse screening length in this channel) 
is about ~ 2(7^T. This indicates that non-perturbative effects are rather large and 
that presumably Cm ~ 2. If true, a perturbative analysis of the effective potential 
in a parameter range relevant for cosmology has low chances to remain if force after 
non-perturbative effects are taken into account. 

In order to clarify this question, one has to use non-perturbative methods for the 
study of the effective potential. They are provided by lattice Monte-Carlo simulations. 
Some preliminary data derived with the use of 4-dimensional lattices is already available 
T^-|l^. There is an indication fl^ that the strength of the first order phase transition 



is enhanced in comparison with perturbative calculations. 
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The purpose of this paper is to combine perturbative and non-perturbative methods 
for the study of the finite T electroweak transition. Namely, with the use of dimensional 
reduction ||2y]-[^ one can first derive a 3-dimensional effective action with perturba- 



tively calculable coefficients, summing up thus the effects of the Debye screening which 
are well understood. Then, one can simulate the phase transition with lattice Monte- 
Carlo methods. 

There are a number of advantages of this method in comparison with 4-dimensional 
lattice simulations. First, it separates the physics in which we are reasonably confident 
(Debye screening) from the unknown non-perturbative 3-dimensional physics. In other 
words, any 4-dimensional simulations contain perturbative noise which has nothing to 
do with non-perturbative effects determining the order of the phase transition. More- 
over, perturbation theory signals that the first order nature of the electroweak phase 
transition is an exclusively 3-dimensional phenomenon, since the so-called cj)^ term in 
the effective potential (this term induces a jump of the order parameter) arises due to 
infrared singularities in loop integrations at high temperatures. The last, but not least 
advantage is that a lattice simulation of 3-dimensional theories is less time consuming 
and more transparent from the point of view of scaling behaviour. 

The paper is organized as follows. In Section 2 we derive the 3-dimensional effective 
action for our system. In Section 3 we compute the effective potential in the continuum 
3-dimensional theory and show how the first order nature of the phase transition arises 
in the loop expansion of the effective potential. In Section 4 we study the effective 
potential for the lattice version of the theory which is then compared in Section 5 with 
Monte Carlo simulations. Section 6 contains our conclusions. 



2 The effective action 

Finite T field theoretic systems are characterised by fields defined over the interval 
0<r</9 = l/Tin imaginary time and extended beyond this region by the con- 
dition of periodicity (antiperiodicity) with the period (3 for bosonic (fermionic) fields. 
If the action is expressed in terms of the Fourier components, the quadratic terms 
are of the type [(2'7rnT)^ -|- k^]|y4(n, k)p, where y4(n, k) is a generic bosonic field and 
n = —oo, . . . , +00. At high T and k < 2tiT the nonstatic modes A{n ^ 0, k) are 
thus suppressed by the factor {2nnT)^ relative to the static y4(0,k) modes. The idea 
of dimensional reduction is to exploit this suppression by integrating over the non- 
static modes in S[A{n = 0, k), A(n ^ 0, k)] and deriving in this way an effective action 
5'efr[A(0, k)] for the dominant static modes. We shall below carry this out for the four- 
dimensional SU(2) + a fundamental Higgs theory, taking g' = in the electroweak 
sector of the standard model. The effective theory then is a three-dimensional SU(2) 
+ a fundamental Higgs + an adjoint Higgs model with well defined T-dependent coef- 
ficients. 

For fermionic fields the square of the inverse propagator is [(2n+ l)7rT]^ + k^ and all 
modes are suppressed at large T. However, the fermionic fields will enter by changing 
the coefficients of the effective action of the bosonic static modes. Their effect can thus 
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also be studied in this framework. However, as our calculations so far are only on the 
one loop level, we can only include a top quark of mass less than 79 GeV, given by one 
loop stability [p5[] . 

The starting point is the action of the 4d SU(2) + fundamental Higgs model 

^[A;:(r,x),0,,(r,x)] = J^r J ^/Mi^^ + P^0)^Pm</>) + /^VV + A(0V)'}, (1) 

in standard notation. Possible fermionic terms are not shown explicitly. 

Including 1-loop corrections and terms not damped by powers of 1/T the dimension- 
ally reduced effective action becomes 

^efrlA^x), AS(x),0,(x)] = i| d^^Ft^Ft^ + i(AAo)'^(A^)'^ + (A0)nA0) + 



1 

+ 2 



+ ^(1 + ^ - + 



.J(l + i + ^)/T2-(4 + l)/TS.]AgAS + -^' - ' ^ 
2 ^ ^8^ 16^ 2 8m^ ' 

1 



127r2^ 16 

A(0t0)2 + 



+ 4^7%"^o</'V^ (2) 



where Sc is the integral 



(27r)3p2 

depending linearly on the cutoff. In perturbation theory it cancels against 1-loop diver- 
gences, as shown explicitly below. If one wants to study the system nonperturbatively 



with a finite cutoff it must be included for correct continuum limit [^. We shall fix 
the parameters of the effective action by taking g = 2mwiV^GF)~^^'^ = 2/3, rriiy = 
80.6 GeV and giving the value of mn- Then X/g"^ = m|^/(8m^). 

The effective theory thus is a 3d SU(2) + fundamental Higgs + adjoint Higgs theory 
with coefficients depending on T, Np, the coupling constants and the cutoff. All the 
three kinetic terms, the two T = potential terms for and the last Aq — (p coupling 
term arise from eq. (|l|) by naive dimensional reduction (taking fields constant in r). The 
general structure of the 1-loop quadratic potential terms for Aq and is CiT^ — C2TX 
cutoff, where the first term is the usual 4d 1-loop screening mass and the second term 
arises from the exclusion of the n = term in the 1-loop integral (since one must only 
integrate over the nonstatic modes). Equivalently, this term is the mass counterterm 
of the superrenormalisable 3d theory. The contributions of the various fields to these 
terms in eq. (0) are ordered so that in the coefficient of 0V fii'st come the terms with 
Af in the loop, then Aq, then and last the fermion loop (only in the terms). In 
the quadratic and quartic Aq terms the ordering is A^, 0, fermions in the loop. The Aq 
loop gives a contribution of order g^ and is neglected here. Also neglected are small 
1-loop corrections to the quartic term. 

In addition to the mass counterterm eq. (§) there in 3d actually also is |2^ a loga- 



rithmic 2-loop counterterm ~ A^T^/(167r^) / (ip/|p|. This is numerically negligible in 
our case. 
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More generally, eq. @ contains terms with higher powers of fields and their deriva- 
tives, allowed by the gauge (BRS) invariance of the 3d theory. For example, the term 
q^^AqAqY /T"^ would appear. If we take the correlation length of the Aq field to be 
l/{gT) we can estimate from the quadratic term that < A'^Aq >~ gT"^. Thus the ratio 
of terms with consecutive powers of A^Aq is 0{g^) <^ 1. 

3 The one-loop effective potential in continuum. 

To study the theory defined by the action eq. (||) in perturbation theory, we compute 
the 1-loop correction Viioop(0, ^o) to the tree potential defined by eq. (|^). One finds 
before regularisation that 

VnooM,Ao)=Tj ^{21og(p^ + /A^ + l/02)+log(p2 + l^V) + 

+^ log(p2 + /i^ + X<P' + ^g'Al) + log(p2 + ml + ^g'cp' + X^A^) + (4) 

+^ log{p^ + p2[/i2 + ml+ (3A + + i^g' + 3Xa)AI] + 

+ 3A0^ + \g'Al){ml + 1/02 + 3XaAI) - \g'Al<p'}], 

where 0^ = 2(f)'' (p and A^ = AqA^. The first two terms come from the Ai loop and the 
remaining ones from the coupled Aq, loops. Using 

/ ^^'^ '"'^'^^^ ^ '"'^'^^^^ " ^f'"' + 

one finds that the Sc terms in eq. (|^) and eq. @ cancel (taking into account the Aq 
loop term —^5XaTT,cAI neglected in eq. (0) for smallness). The finite terms give the 
quantum correction to the tree potential. The resulting 1-loop improved potential is 
then 

V{<P, Aq) = 1/.202 + ^mlAl + \x{<pr + \UAir + l9'Al<P' - 

-£{2/(l02 + A^p + i,V+ (6) 
+ {ml + \g'<P' + XaAD'/' + l{ml + \g^' + 3A^A^)3/2 + 
+ A0^ + -/Alf' + ^(/i^ + 3A02 + \g'Alf^], 

where 

We have, for completeness included the small A^ term, which will again be neglected 
from now on. 
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Studying the minima of V{(f), Aq) one sees that they can be driven to nonzero values 
of Aq only if the negative term is large. This demands that at least g > n^S/S which 
is beyond the domain of validity of this calculation. We conclude that the minimum is 
always at Aq = 0; no condensate is formed. Basically this is due to the large value of 
rriD and the small value of the correction, ~ — T/Gvr. 

For Aq = and \ <^ the corrections to the effective potential have the form 

It is the first term which gives rise to the first order phase transition in perturbation 
theory. Note that the effective 3-dimensional theory correctly takes into account Debye 
screening, which decreases the magnitude of the cubic term relative to the naive loop 
expansion by the factor 2/3 |0. The second term corrects the 7 in eq. (|^ by the term 
(— 3/167r)^2/3(yf^. With this input the equation 

gives the perturbative result for the transition temperature Tc and the lower end Tq < Tc 
of the metastability range (the high T phase does not exist for T < Tq). Numerical 
values for these as well as for some other relevant quantities as calculated from 1-loop 
perturbation theory [|^, [|IU| are given in Table 1. For g = 2/3 the relevant formulas 
are ^{T,)T, = 9TtmH/mw, Tjmw{T,) = 9nml/{2ml,), a/Tj^ = 2ml^,/ {24371^ m%), 
L/T^ = 2m^[l/6 + mj,/{18ml.) - ml,/ il8n^ml)]/ (gn^mj,). 



rriH 


To Te r+ 


m) 




a 


L 


35 


90.89 95.24 95.83 


12.28/Tc 


2.67/Tc 


0.017T3 


0.085T^ 


80 


182.34 183.55 183.71 


28.1/T, 


13.9/T, 


0.00028T3 


0.0044T^'' 



Table 1: Values of Tc, the lower (Tq) and upper (T+) ends of the metastability range, 
the correlation length ^ for the Higgs field at T = Tc, the gauge field correlation 
length l/mw{Tc), the interface tension a and the latent heat L as calculated from 
1-loop perturbation theory for mn = 35 and 80 GeV, for g = 2/3 and for no fermions 
included. 



4 The lattice action and effective potential. 
4.1 Lattice action. 

To latticize eq. (^ we go over to the matrix representation Aq = A^T"-, T"- = ^aa, 
— >• $ = (00 + iai(j)i)/\/2 and rescale the fields by 

igaAQ-.AQ, $ ^ (10) 

V a 2 

5 



where a is the lattice spacing. The lattice action on an A^^ lattice then becomes 
^ = /?GEE(l-^TrP,,) + 

X i<j ^ 

+ ^/?GEE[Tr^o(x)f/r^(x)Ao(x + Of/,(x) -TrAg(x)] + 

X i 



+ ^|lOS(iV3 



b + Np 4 1 1 



G 



2TrA^(x) + 



g'^G,. 1 



^ 37r2 ^ 16 8 '^2 °^ " 

+I3h E E[^Tr $t(x)$(x) - ^Tr <l>t(x)?7,(x)<l>(x + i)] + 

+ - 2/?« - 3/5H)^Tr $t(x)$(x) + /5«(lTr $t(x)$(x))^] + 

-^/?^^E[^TrA^(x)iTr<|.t(x)<|.(x)], 



where 



L-l 

E 



(12) 



4A^3 L-^^ sin^(7rni/iV) + sin^ (vmz/iV) + ^vc^iT^n^jN) 

(the term with all = is omitted; numerically Y.{N^) = 0.224605, 0.233942, 0.238633, 
0.243326, 0.252731 for = 8,12,16,24,oo, respectively) and 



'G 



4 1 
5(2 Ta 



i5n - -XTaP, - 



(13) 



AT) 



2(1 - 2f3R - 3f3H] 



In eq. (|T3|) yU^(T) is the coefficient of the (f^cf) term in eq. (|^) with the continuum Sc 
replaced by S(A^'^)/a. Introducing this relates [3h and T for given values of g, mw, i^h 
(and possibly of mtop): 



m 



H 



4r2 



4 



G 



3--^+2 



f3H 9 



/3h 8m2j,/3G 2(3g' 3m 



|-)E(iV3 



2 

top 



2 V16 ■ 2 8m^ 8m^/' 



(14) 

The curve T = T{(3h) is plotted in Fig.l for a few parameter values. The curves 
shift to the right (left) with increasing mn and S(A^'^) (/3g)- For example, the values 
of Ph corresponding to T = Tq and T = oo are 



/3h(To) = - + (1 + 



H 



3m 



S(iV^) 
/5c 



f0(A 



W PG 



G )^ 



m 



H 



(15) 
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Eq.(0) takes correctly into account the divergent mass counterterms of the 3d theory. 
In addition to these, there will be finite renormalisations of the coupling constants g 
and A. Due to infrared divergences getting more and more serious at higher orders 
T7[] these renormalisations - which will depend on the lattice size N - cannot be 



calculated perturbatively. We shall observe that the numerical results will follow the 
"constant physics" curve eq. (|T^, with dependent mass renormalisation but N 
independent bare values of g and A, rather well. However, some more dependence, 
clearly attributable to finite renormalisations of the coupling constants, will remain. 
In any case, all finite size effects will be very difficult to control. For example, it is 
not known how to include the constant mode (the nj = term in eq. (|T^) ) term in 
perturbative calculations. 



4.2 Effective potential on the lattice. 



Monte-Carlo lattice simulations should be compared with the perturbative calculations 
of the effective potential on the lattice rather than with continuum expressions. To get 
a lattice generalization for the effective potential one can just change the integration 
over momenta in eq. (^ to a finite sum over the discrete momenta Pi = {2tt / aN)ni, rii = 
0, N — 1. Including only the Ai loop term, which is relevant for the first order nature 
of the transition, the lattice effective potential becomes 



latt 



where 



^7(T^ 



d 



3T 



E 

ni=0 



log 1 + 



(W/4)^ (gaH^)' 



d 



d 



sm 



(Tini/N) + sin^(7rn2/A^) + sin^(7m3/A^). 



, (16) 



(17) 



Note that in this sum the term with rii = must not be included not only since we 
cannot handle it but also by definition since the effective potential is constructed for 
x-independent and thus constant configurations must not be integrated over. 
To study the order of the transition it is convenient to study the zeroes of 



dVv 



latt 



N-1 



7(T2-T2) + A0^ 



3/ ,2 ^ 1 



N3 



d[d + {ga(j)/Ay 



In the continuum limit a — 0, — > cxd, aN — >■ oo the last term becomes — (3/16'7r)5f^T0. 
For T = Tq one then is solving A0^ — {3/16'n')g^T(j) = 0, which trivially leads to 
a second minimum. This also exists for some temperatures above Tq, up to T = 
T+. The case of finite lattices is completely different. It is clear from the lattice 
expression for the effective potential that it is an analytic function of 0^ at least for 
0^ < (IG/^f^a^) sin^(7r/A^) and , therefore, no term ~ —0 can appear. Nevertheless, 
a second minimum can exist for sufficiently large lattice sizes. The condition for this 
simply is that the equation 



A 



128 



N-1 



E 



'q d[d + {ga(j)/Ay 



(19) 
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have a solution. The right hand side decreases monotonically when increases and we 
can take = 0; if a solution appears it appears first at = 0. Inserting aT = 4/ (g'^Pc) 
and A = g'^m'jj / {8m^r) gives the relation 

The approximation is the result of an explicit numerical calculation of the sum. For 
the value (3c = 20 used in our numerical simulations one sees that for m^^ = 35 GeV 
a first order transition appears for > 29 and for rriH = 80 GeV for > 160. The 
use of so large lattice volumes is not possible in practice, and we at this stage confined 
ourselves to smaller N, N < 20. If perturbation theory works well, one should not find 
any signal of a first order phase transition in the Monte-Carlo simulations with these 
small lattices. 

It may be of interest to note that in the limit of large lattices (a = constant, —>■ oo) 
the last term in eq. ([iBD can be written in the form ||29l 



Jo 



^{l-exp[-{ga^/4fy])-{ga<P/^f 



exp{-3y)l',iy), (21) 



with the aid of the Bessel function Iq. Numerically this is very close to the continuum 
result, in particular, the 0^0 limits are the same. In other words, the critical 
temperature as well as the metastability range given in Table 1 are practically the 
same for very large N and for the continuum. 



5 Lattice simulations 

Our choice of parameters for the simulations is motivated as follows. We have already 
fixed that g = 2/3 and mw = 80.6 GeV. Since our aim is to test the validity of 
perturbation theory we choose ttlh = 35 GeV, which makes A small = 0.0105 but is not 
too close to the vacuum stability limits of somewhat less than 10 GeV. For comparison 
we also choose ttih = 80 GeV (A=0.0547). We also do not expect fermions (except 
perhaps for the top quark) to qualitatively change the nature of the transition and thus 
choose Np = 0. Light fermions could simply be included by changing the numerical 
value of the coefficients as in eq. (|^). Finally, we choose Pg = 20, which basically fixes 
the lattice spacing a in physical units via the first equation in eq. (p!B|): a = 0.45/T. 



As should, this is smaller than the thermal distance scale 1/T, the average distance 
between particles. Equivalently, in order to describe correctly effects associated with 
magnetic sector of the 4-dimensional theory one must have a -C I/ttim- With the 
use of Cm ~ 2 and Pg = 20 we get a ^ 0.5/mM- For this a the perturbative Higgs 
field correlation lengths in Table 1 are 27a {mfj = 35) and 42a [rriH = 80) while the 
gauge field correlation lengths are 6a (mi^=35) and 14a {rriH = 80). These are so large 
that it is realistically not possible to fit an interface between two phases (broken and 
unbroken) in the lattice. 



8 



5.1 The update algorithm. 



Because the lattice system described by the action eq. (0) has several qualitatively 
different components, we used a wide mixture of update algorithms to obtain good 
performance. The SU(2) gauge field was updated as follows: first, we combined the 
plaquette action (first term in eq. (0)) and the $-field hopping term (fifth term) 
to form a local SU(2) action of the form /3'j;^{x)TT Xi{x)Ui{x) , X G SU(2). Using 
this action, new link matrices were generated with the Kennedy-Pendleton heat bath 



method The adjoint field hopping action 5'hopp(f^, ^o) (second term), which is 

quadratic in Ui, was then taken into account by accepting or rejecting the new link 
matrices with the Metropohs method - that is, with the condition exp[S'hopp(f^°'^'*, Aq) — 
5'hopp(^°^'^; ^o)] > where r is a random number from an uniform distribution between 
and 1. The acceptance rate for the Kennedy-Pendleton heat bath was ~ 99.5% and 
for the accept/reject step ~ 95%. 

The length of the adjoint field Ra = {AqAqY^"^ was updated with the Metropolis 
method, while the colour space direction, which appears only in the hopping term, was 
updated with an SO (3) heat bath. Similarly, the fundamental Higgs field was divided 
into radial and SU(2) parts: $ = RV, R>0,V e SU(2). The l^-field, which appears 
only in the hopping term (the fifth term in eq. (|TTp), was updated with the Kennedy- 
Pendleton heat bath algorithm, and the length R was updated with the Metropolis 
algorithm. 

The evolution of the gauge and Ao-fields in the simulation time was very rapid 
compared to the evolution of the $-field. It is crucial to make the Higgs field update 
as effective as possible. Because of the large Pg the gauge background of the Higgs 
field is very flat, and it is plausible that one could construct an effective multigrid or 
cluster update algorithm. 



5.2 The results. 

Results of the simulations are shown in Figs. 2-6. In Fig. 2 we present the probability 
distribution of the order parameter L = Tr y^(x)f/j(x)l^(x + i) for niH = 35 GeV and 
different lattice sizes. The two peak structure characteristic of a first order phase tran- 
sition is clearly seen on these 8^ and 20^ lattices, as well as on intermediate 12^ and 16^ 
lattices (not shown). The change of the curves with the lattice volume is qualitatively 
consistent with what one would expect from a first order transition, namely the two 
peak structure is more distinguished on the large lattices. Note that the peak positions 
do not depend on the lattice size which indicates that the finite volume effects are not 
substantial. 

The first order nature of the phase transition can be also seen in Fig. 3 where the 
simulation time evolution of the order parameter L is shown. Here initially the system 
was confined in the unbroken phase with small value of L, then it jumps to the phase 
with broken symmetry and stays there quite a long time. These jumps then continue 
in Monte Carlo "time". 

An important criterion for the order of the phase transition is the finite size scaling 
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of the second moment of the order parameter L [^. This is shown on Fig. 4 for 
various lattice sizes. The continuous curves were obtained by combining the individual 
runs (performed at various (3hs) with the multiple histogram (Ferrenberg-Swendsen) 
method We observe that the second moment of L, as a function of [3h, develops 
a narrow peak as the volume is increased. The height and the location of the peak 
have well-defined infinite volume limits; this behaviour is characteristic for systems 
exhibiting a first order phase transition. 

Analogous data for a more heavy Higgs {niH = 80 GeV) is presented on Figs. 5-6. 
In this case the correlation lengths for the Higgs and W bosons are larger than the 
lattice size and strong volume dependence is observed (Fig. 6). No definite conclusion 
can be made for this case with the present lattice sizes. 

According to Fig. 2 the transition for niH = 35 GeV takes place at about Ph = 
0.34010 roughly independent of the lattice size. We thus do not observe precise scaling 
according to eq. (p!4D, which would demand that the value corresponding to some fixed 
temperature change according to eq. (|I^) with lattice size A^. We ascribe this to the 
nonperturbative renormalisations of g and A discussed above. If we extrapolate to the 
curve corresponding to A^ = cxd, we see that [3h = 0.34010 corresponds to Tc ~ 85 GeV, 
somewhat but not much below the continuum perturbative value of Tc = 95 GeV. For 
ruH = 80 GeV the results are less conclusive, but for the largest lattice studied the 
transition takes place at (3h = 0.3418. The uih dependence of eq. (|15|) is rather well 
reproduced, but so far it is impossible to make a definite conclusion concerning A^ 
dependence. Anyway it is suggested by Fig.l that again the observed is less than 
the perturbative value of 184 GeV. 



6 Conclusions 

We have performed a combined analytical and numerical study of the finite T elec- 
troweak phase transition. First those degrees of freedom which - with reasonable 
degree of confidence - can be treated perturbatively were analytically integrated over 
and an effective action in the remaining degrees of freedom was derived. This effective 
theory is a T = 3d SU(2) gauge field -|- adjoint Higgs + fundamental Higgs theory 
with known coefficients depending on T and, quite essentially, on the cutoff of the 
theory. This bare effective theory was then latticized with 1/a, a = lattice spacing, as 
the cutoff and the nonperturbative degrees of freedom were treated numerically with 
Monte Carlo techniques. 

In this first numerical application our aim was to study the structure of the theory 
and a fairly small Higgs mass, rriH = 35 GeV was mainly discussed, with some results 
given also for tuh = 80 GeV. A first order transition was clearly seen at least for the 
smaller Higgs mass and the numerical value of is rather close to (but less than) 
the value obtained from 1-loop perturbation theory. For this numerical agreement the 
computable cutoff dependence of the bare effective theory was crucial. With increasing 
rriH the various characteristic lengths in the problem increase rapidly and the problem 
becomes numerically more and more difficult. 
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Our lattice Monte Carlo results thus clearly indicate that there is a first order phase 
transition, at least for small values of the Higgs mass. A very intriguing part of the 
result is that we see a first order phase transition also on sma// lattices where, according 
to the perturbative calculations in Section 4.2, the phase transition must be of the 
second order. A number of tests including the lattice volume dependence of the order 
parameter, temperature metastability range, etc. indicate that the first order character 
of the phase transition is not a lattice artifact. These results indicate that so-called 
0^ term is not the only source for the first order character of the electroweak phase 
transition and that non-perturbative effects are important as well. The physical nature 
of these non-perturbative effects should be presumably related to the confinement in 
3-dimensional gauge theory. Unfortunately, the complete solution of the problem of 
the electroweak phase transition by non-perturbative lattice methods requires huge 
lattices. 
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Figure 1: The relation between T and (3h as given by eq. (|T^) for ttlh = 35 and 80 
GeV, Pg = 20 and Np = 0. The arrows denote the perturbative values of as given 
in Table 1. 
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Figure 4: The dependence of the second moment of the order parameter L on for 
different lattice sizes and uih — 35 GeV. 
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Figure 5: The distribution of L = Tryt(x)C/i(x)y(x + i) for tuh = 80 GeV and a 16=^ 
lattice for different values of (3h- 
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Figure 6: The dependence of the second moment of the order parameter L on Ph for 
different lattice sizes and niH — 80 GeV. 
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